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ABSTRACT 

Chemical models used to study the chemical composition of the gas and the ices in the 
interstellar medium are based on a network of chemical reactions and associated rate coefficients. 
These reactions and rate coefficients are partially compiled from data in the literature, when 
available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public 
gas-phase network first released in 2012. In addition to a description of the many specific updates, 
we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as 
compared with the results of the previous version of our network, kida.uva.2011. 

Subject headings: astrochemistry astronomical databases: miscellaneous ISM: abundances ISM: 
molecules 
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1. Introduction 

Astrochemical models have been developed 
over the years to study the chemical compo¬ 
sition of the interstellar medium (ISM) (see 
Agundez fc Wakelaml2013i and references therein). 
Current models follow the abundance of hun¬ 
dreds of specie s thro ugh thousands of reactions 
( McElrov et al.1 l2013f) . Although the processes 
at the surface of the grains, such as sticking, 
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surface reactions, an d desorption, are crucial 
Ij Herbst fe Yates! 12013 ), gas-phase reactions are of 
prime importance and poor estimates of rate coef¬ 
ficients can induce large errors in th e predicted 
abundances, as occurs for OCS ( |Loison et al 
20121) . The rate coefficients for chemical reac¬ 


tions under the ISM conditions of low density and 
temperatures from about 10 K to hundreds of K 
can be studied in the laboratory for som e systems 
or studied theoretic ally for others (Smith 12011 : 
Larsson et ah ia. Because each study on a 


specific reaction can take from a few months to 
several years, such studies are best undertaken 
for systems identified as the most important ones 
( Wakelam et al.[201(il . For other systems, general 
rules can be used (with caution) to estimate the 
rate coeffi cients by similariti es to other studied 
reactions ( Loison et al • 2014dh Even for systems 
studied in detail, it is uncommon to explore tem¬ 
peratures as low as 10 K, so that some estimate or 
extrapolation of the temperature dependence may 
be necessary but must then be performed with 
caution. 

The KIDA online database!]] is a compilation of 
rate coefficients coming from various sources in¬ 
cluding published papers, and other databases or 
data sheets, which are of interest for the ISM or 
the modeling of planetary atmospheres. For “im¬ 
portant” reactions, the appropriate KIDA experts 
can be asked to provide data sheets with recom¬ 
mendations on the rate coefficients to be usec0. 
KIDA contains about 70 such data sheets and the 
number is increasing with time. The list of ex¬ 
perts in physical chemistry is provided in Table [1] 
with their current coordinates. To help astro- 
chemical users, we also provide a subset of chemi¬ 
cal reactions, extracted from KIDA, to be used to 
model the chemistry in the ISM. The first release 
of this network , kida .uva.2011, was reported in 
Wakelam et al.l ( 2012 ') (hereinafter W12). We are 


now presenting in this paper an updated version, 
labelled kida.uva.2014. The chemical updates of 
the network are described in Section 2 while the 
impact of the new chemistry is discussed in Section 
3. 


1 http: //kida.obs.u-bordeauxl.fr/ 

2 http: //kida.obs.u-bordeauxl.fr/datasheets 


2. Chemical updates 

This new network is based on the previously 
published kida.uva.2011 network (W12) with the 
updates posted in the KIDA database between Oc¬ 
tober 2011 and October 2014. The main updates 
of the network come from: 


- Reactio ns from the high te mperature net¬ 


work of Harada et al. ( 2010h , 


Photo rates from 

van Dishoeck (1988); 

Roberge et al. (1991 

): van Dishoeck et al. 

(2006h; van Hemert & van Dishoeck ( 

2008), 


- Rate coefficients for the OCS chemistry from 


Loiso n et al.1 (2012), 


- New data sheets (IWakelam et al.ll2013l and 
many others), 


- A review of the HC N/HNC chemistry by 
Loison et al. ( 2014bl h 


- A revie w of the carbon chemistry bv lLoison et al 
(l2014cJ) . 


- Branching ratios for reactions forming C^y^_ 10 , 
CW - 4 H(°’+), and C 3 h'°’ +) from lChabot, et al. 

E3). 


In addition to these updates, a number of individ¬ 
ual reactions have been added or updated from a 
list of studies given in the Appendix. New species 
have been added to the network, and are listed 
in Table [2] A total of 446 rate coefficients have 
been changed and 1038 new reactions have been 
added. The final network contains 489 species 
composed of 13 elements and 7509 reactions. In 
addition to the gaseous species, neutral grains, 
negatively charged grains, grain-surface electrons 
and atomic hydrogen are included. The types of 
reactions and the formulas used to parameterize 
the rate coefficients are the same as described 
in W12. The complete network is available at 
http://kida.obs.u-bordeauxl.fr/uploads/models/kida.uva.2014.zip 
with a list of species and the bibliographic refer¬ 
ence in bibtex format. 


As for the previous version, kida.uva.2011, we 
expect the network to be used for the chemical 
modeling of the interstellar medium at temper¬ 
atures between 10 and 300 K. Despite the fact 


that we added reactions from Harada et al. ( 2010l ) 
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for higher temperature conditions, including re¬ 
actions with activation energy, and reverse reac¬ 
tions, we have not done any additional checks on 
its appropriateness for temperatures larger than 
300 K. For some of the reactions, the parameter 
7 used in the temperature-dependent Arrhenius- 
Kooij expression for the rate coefficient - k(T) = 
a(T/300) /3 exp(—y/T) - is negative. This sign 
change will produce strongly incorrect results if 
the calculation of the rate coefficient is extrap¬ 
olated outside the temperature range of validity. 
Those reactions are listed in Table [3] Not extrap¬ 
olating outside of the range of validity can also 
have consequences. In fact, the temperature de¬ 
pendence of rate coefficients is usually valid over a 
certain range of temperature, which can be a cru¬ 
cial piece of information. However, in many cases, 
the information is only partial and could lead to 
large errors if not carefully considered. This is par¬ 
ticularly true for H-atom transfer reactions where 
quantum mechanical tunnelling results in strong 
departure from Arrhenius be havior at low tem¬ 
peratures ( Tizniti et a, 1.1 EffiTlt ). As an example, 
consider the reaction N + H 2 —> H + NH, which 
has a very large activation energy. This reaction 
was already present in the previous OSU database 
( Harada et alJl2010r i without any indication of the 
temperature range and, it is likely that common 
users would extrapolate the rate coefficient down 
to 10 K. In KIDA, we incl uded the original rate 
coefficient from Davidson & Hanson ( 19901 1 with 
its measured temperature range of 1950 - 2850 K. 
The rate coefficient computed at the minimal tem¬ 
perature (1950 K) would then be 4 x 10 ” 13 cm 3 
s ” 1 (much larger than the 10 K extrapolation of 
4 x 10 ” 82 cm 3 s” 1 ). If only the measured temper¬ 
ature range for N + H 2 (1950 - 2850 K ) were in¬ 
cluded in the network, the reaction rate coefficient 
for temperatures lower than the lowest measured 
temperature (1950 K) would be programmed to 
remain at the measured value for 1950 K. Like¬ 
wise, the rate coefficient at temperatures higher 
than 2850 K would be programmed to remain at 
the 2850 K value. In dense cloud conditions where 
both N and H 2 are abundant, even the low 1950 
K rate coefficient would be large enough to have 
a significant effect on the chemistry while the rate 
coefficient of this reaction at 10 K is clearly neg¬ 
ligible. For this reason, in kida.uva.2014, we have 
used the original rate coefficient and reference but 


we have changed the temperature range to extend 
to 10 K to allow extrapolation. For reactions be¬ 
tween ions and non-polar neutrals with measure¬ 
ments only at 300 K, we have assumed that the 
rate coefficients are independent of temperature, 
as in the Langevin approximation, and have used 
the measured total rate coefficients and product 
branching fractions over the temperature range 
10-300 K. 

3. Modeling results of the new network 
3.1. Comparison with kida.uva.2011 

To quantify the impact of the updates done on 
kida.uva. 2011 , we have compared the abundances 
computed with the chemical code Nahoon using 
kida.uva.2011 (see W12) and kida.uva.2014 with 
physical parameters typical for dense clouds: a 
temperature of 10 K, a total “proton” density, 
defined by the relation «h = n(H) + 2n(H2), of 
2 x 10 4 cm” 3 , a cosnric-ray ionisation rate for H 2 of 
10 ” 17 s ” 1 and a visual extinction of 30. The gas- 
phase abundances of the species are computed as a 
function of time starting from an atomic composi¬ 
tion except for hydrogen, which is initially molec¬ 
ular. The elem ental composit i on ut ilised is the 
same as used by iHincelin et al' ( 201llh except for 
elemental oxygen, which has an abundance with 
respect to H of 2.4x 10” 4 . Tables[4]and Ogive lists 
of abundant species, defined as those with abun¬ 
dances computed with at least one network larger 
than 10” 12 compared with the total proton den¬ 
sity. The lists contain only those abundant species 
with abundances modified by more than a factor of 
2 at 10 5 and 10 6 yr with the change in reaction net¬ 
work. For each of these species, we also give the ra¬ 
tio of the abundance computed with kida.uva.2014 
to the one computed with kida.uva .2011 in loga¬ 
rithmic terms. In Fig. [lj we show the abundances 
of a selection of species as functions of time com¬ 
puted with the two networks. 

At 10 5 yr, the species that show a difference 
larger than two orders of magnitude are HC 2 N + , 
HC n N (with n=4, 6 and 8 ), CH 2 NH and CCS, 
and they are all decreased by the new network. 
The HC 2 N + abundance is decreased by the newly 
added reaction H 2 + HC 2 N + —> H + CH 2 CN + 
while the HC n N species and CH 2 NH are decreased 
by new d estruction reactions w ith atomic carbon, 
all from Loison et ah ([2014b). The strong de- 
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crease of the CCS abundance is due to the in¬ 
troduction of the destruction reaction O + CCS 
—> CO + CS, which possesses a rate coefficient 
of 1CT 10 cm 3 s -1 independent of temperature 


( Loison et al. 2012h . The rate coefficient has been 
deduced from the simiiar O + CCO reaction for 
which the rate coefficien t has been measured to 
be 9 x 10 -11 cm 3 s _1 ( Shackleford et all Il972{ ; 


Bauer et al.lll985h at room temperature. The high 
value of the rate coefficient is compatible with the 
fact that the O + CCS reaction is a radical-radical 
reaction which usually shows no (or a submerged) 
b arrier. Note that this reaction was also included 
in lHarada et al . ( 201 01 but with an activation bar¬ 
rier based on the rate coefficient from O + CS. 

A significant fraction of the molecules listed 
in Table 0] are either nitrogen-bearing species or 
carbon-bearing species because a large number of 
the rate coefficients for these species have been 
changed, as mentioned in Section [2] The abun¬ 
dances of large carbon-chains have also been re¬ 
duced. At 10 6 yr, the list of species is much smaller 
than at 10 5 yr and the differences are smaller 
than at earlier times, partly because the abun¬ 
dance of atomic carbon decreases with increasing 
time. Among the species that present a variation 
of more than a factor of 10, only NO was not 
already present in Table 3. This molecule is in¬ 
creased by approximately a factor of 10 because 
of the decrease in the N + NO rate coefficient as 
recommended by the KIDA experts in a datasheet. 

3.2. Comparison of model predictions 
with observations 

Although this is not the primary goal of this 
paper, we have compared the model results, dis¬ 
cussed in the previous section, with th e observa¬ 
tions in TMC-1, as listed in Agundez &; Wakelaml 
( 20131 . Table 4). We initially used only gas-pliase 
chemistry and in a second attempt, we considered 
grain-surface chemistry as well. For this second 
test, we used the Nautilus gas-grain code as de¬ 


scrib e d in many previous p apers (lSemenov_etal 
201fll : R.eboussin et al. 20l4 Loison et al.1 2014cl^ 


The grain surface network was adapted to in¬ 
clude the new species introduced in kida.u va.2 014 
but is basically the same as in ILoison et al 


( 2014ch . It consists mainly of diffusive Langmuir- 
Hinslrelwood processes with adsorption and des¬ 
orption mechanisms. In addition to thermal des¬ 


orption, Nautilus includes cosm i c-ray ind uced 
desorption, as described in Hasegawa fe Herbst 
and the desorption by exotherm ic ity o f 


(199- 


surface reactions following iGarrod et al. (1200711 


with the parameter a set to 0 . 01 . 

The physical conditions and initial abundances 
are the same as described in the previous section. 
To compare with the observations, we have sim¬ 
ply computed the differences between the loga¬ 
rithms of observed and modeled abundances us¬ 
ing the two networks, and determined the time 
for which t he mean logarithmi c difference is at its 
minimum ( Loison et alJl2014c l. This logarithmic 
parameter is defined by the sum over all species 
of |log(X mo d/X 0 bs)| /N 0 b s , where X nKH | and X Q bs 
are the modeled and observed abundances while 
N 0 bs is the number of observed species (57) in 
TMC-1 used for the comparison. Considering only 
gas-phase chemistry, the general agreement with 
the model is slightly worse with the new network. 
When surface chemistry is considered, we obtain 
a similar degree of agreement but for different 
times. We now only discuss the comparison us¬ 
ing the full gas-grain networks. Fig. [2] shows the 
mean logarithmic difference as a function of time 
for both networks. We have ignored the species for 
which only upper limits have been constrained. At 
the best times, for both models, the mean differ¬ 
ence between modeled and observed abundances 
is smaller than a fac tor of 10. As already noted by 


Loison et al, ( 2014cT) . we now obtain a reasonable 


agreement with the observations in TMC-1 using 
a C/O elemental ratio smaller than 1. The best 
time is 2.5 x 10 5 yr using kida.uva.2011 while it is 
approximately 1.3 x 10 6 yr using the new network. 
For these best times, we have plotted in Fig.0]the 
difference between modeled and observed abun¬ 
dances for individual species. The difference be¬ 
tween modeled and observed abundances is de¬ 
creased for about half of the species and increased 
for the other half when comparing the two net¬ 
works. This explains the fact that the general 
agreement is not changed. 

The two species for which the changes are 
the strongest are CH 3 CHO, for which the agree¬ 
ment with the new model is worse, and CCO, 
for which agreement with the new model is bet¬ 
ter. Fig. [4] shows the abundance of these two 
species as a function of time predicted with the 
two networks. The TMC-1 abundance of CCO is 
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3 x 10~ n (compared to the total proton density) 
( Ohishi & ; Kaifu| 1998 ) while tha t of CH 3 CHO is 
3 x 1Q ~ 10 (jMatthews et ah 1985). Both modeled 
abundances are smaller than the observed ones at 
all times. The CCO predicted abundance is much 
smaller at 2.5 x 10 5 yr, the best time using the 
older network, than at 1.3 x 10 6 yr. After 10 6 yr 
however, the CCO abundance is also larger with 
the new network because of the larger abundance 
of atomic carbon since CCO is produced by C 
+ HCO —y CCO + H. The better agreement be¬ 
tween the new model and the observed CCO abun¬ 
dance is then due to the change in the best time 
rather than the chemistry itself. For CH 3 CHO, 
the worse agreement with the new model is also 
due to the change in the best time: at later times, 
the CH 3 CHO abundance is smaller. The decrease 
in the gas-phase abundance of this species is how¬ 
ever much stronger with the new network as can 
be seen in Fig. 0] At 10 6 yr, one of the main for¬ 
mation reactions of CH 3 CHO is O + C 2 H 5 —>• H 
+ CH 3 CHO. The CH 3 CHO abundance is then de¬ 
creased with a decrease in the C 2 H 5 abundance. 
The C 2 H 5 molecule is more strongly destroyed us¬ 
ing the new network because of the newly intro- 
d uced reacti on H + C 2 H 5 —> CH 3 + CH 3 from 
Baulch et ah (119921) . The rate coefficient for this 


reaction is 6 x IQ -11 cm~ d s - 1 and was recom¬ 
mended by 1 Baulch et al. (1992) for temperatures 
between 300 and 2000 K. Here we use the same 
rate coefficient down to 10 K. 


4. Conclusion 

The Kinetic Database for Astrochemistry 
(KIDA) aims at providing to the community gas- 
phase kinetic data to study the chemistry under 
the extreme conditions of the interstellar medium. 
From this database, a subset of reactions is ex¬ 
tracted to be directly used in chemical models. 
The first release of this subset, kida.uva.2011, 
was included in W12. After a large number 
of updates, we are now presenting the new ver¬ 
sion: kida.uva.2014. A large number of gas-phase 
species are affected by these modifications at typ¬ 
ical cloud chemical ages (10 5 — 10 6 yr). The mean 
ratio between modeled abundances and the ones 
observed in the TMC-1 dark cloud is a factor of 6 
at best, using a full gas-grain model. This general 
agreement is not much affected by the new gas- 
phase network but the ’’best time” is older than 


obtained with the previous one. 
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S.D.L.P. is grateful to the Institut Universitaire 
de France, while E. H. thanks the NSF and NASA 
for support. The work at Argonne is supported 
by the U.S. Department of Energy, Office of Sci¬ 
ence, Office of Basic Energy Sciences, Division of 
Chemical Sciences, Geosciences, and Biosciences 
under Contract No. DE-AC02-06CH11357. 


5. Appendix 

5.1. References for updates 

The references u s ed to update the network are: 


199^; 

Arthur & Bell ( 

1971); 

Avramenko & Krasnen’kov 

19671); 

Azat van et al. 

1975 I): 

Bauerle et al. 

(1995blal): 


Baulch et al 


Brvukov et al 


( 1992L _ . 

(12001 klChabot et al.l (1201 


I 

al. 


Cohen fe Westberg 


980h : Gustafsson et al.l (2014 ): Maluenc es et al 
( 19931: iMaver et al.l (Il966fh~ Mebel et al. ~ifl996h : 


( 1988 ) : [Mitchell ( 198 4 a] bj) ; Mi yoshi et al.l (119931): 


Murrell fc Rodrigue d 


(20091): Nsruven et. al l (20tk ^Otto^.al 


20121); Daranlot et al. (20131); Diau et al. 

1995 

); 

); 

Dombrowskv & W 

ignerl (1992]); Frank 

1986 

Freund & Palmer 

19771): Fimi et al. (1987) 

Gustafsson et al. 

(2012): Hack et al. 

(19791); Hamberg et al. 

2010 

); 

); 

); 

); 

Hanson & Salimianl (19841): Harada et al. 

2016 

Hemsworth et al. 

19741): Henshaw et al. 

1987 

Herbrechtsmeier 

1973!): Hickson et al. 

2013 


Kim ( 1975lh Kim k. Huntress Jr.l(ll975 1: Leen fc Graff 

1988 ): Lifshitz fe Michael (Il991 ): Llovd (Il974 ): 


oison et al. (20121 2014alb c ^ Macka^.et.jl J 198 


Millar et al. ( 1986 . 19871 1991 ): Miller fe Meliusl 


1986 


Patterson fc Greenel (Il962 ): 


( 1984 . 


19921) 


); iNeufeld fc W olfire 


Payzant et al. 


Peeters fe Mahnenl (Il973l): Roberge et al.1 (19911): 
Rodgers (19961): Schofield! ( 19731): Shannon et al. 


(2014); Singleton fc Cvetanovic 1l988 ): Smith et al 


1992L 119941): Stewart et al. (19891)7 Tsangl 

; Tsang fc HamnsonT i 19861 ) : lTsang fc Herronl 
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Tsuboi fc Hashimotol ^198lh: van Dishoeck 


Bryu kov, M. G., Slagle, I., & R. Knyazev, V. D. 


(1991);. ___ , , „ . . 

(19881: van Dishoeck et, al. ( 20061) : van He mert, fc van Dishoe6kl () 1. J. Phys. Chem. A, 105, 6900 
( 2008^ Vandooren et al.l j 19941): Warnatd (Il984ll: 


Yang et al.l ( 1993l) ~ Zanchet et al. ( 20091 ); Zellner fc Ewigl 


(119881) . 
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Table 2: Species added to the kida.uva.2014 network 


New species 

c 5 o, c 6 n, c 2 h 6 , c 7 o, c 8 n, 

C 9 0, C 10 N, hn 2 o+, hc 5 o+, c 6 n+, hc 6 n+, 
hc 7 o+ c 8 n+, h 2 c 6 n+, c 2 h+, hc 8 n+, hc 9 o+, 
c 10 N+, h 2 c 8 n+, hc 10 n+, h 2 c 10 n+ 


Table 3: Reactions with negative gamma parameters 


Reaction 

T range (K) 

Reference 

F + H 2 —>• H + HF 

10 - 100 

Neufeld et al. (2005) 

OH + OH —»• 0 + H 2 0 

200 - 350 

datasheet 

C + NH 2 -> H + HNC/HNC 

10 - 300 

datasheet 

0 + OH —>■ 0 2 + H 

150 - 500 

datasheet 

H 2 + NH+ -> NH| + H 

10 - 20 

datasheet 

C + OH -> H + CO 

10 - 500 

Zanchet et al. (2009) 

CH + OCS —>■ CO + H + CS 

301 - 667 

Zabarnick et al. (1989) 

N + Si —> SiN + hu 

10 - 100 

Gustafsson et al. (2012) 

H + F -)• HF + his 

100 - 2500 

Gustafsson et al. (2014) 
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Table 4: Species with abundance ratios log(X new /X 0 id) smaller than -0.3 or larger than 0.3 at 10 5 yr, and 
with an abundance in at least one of the models larger than 10' 12 compared with the total proton density. 


Species 

Ratio 

Species 

Ratio 

Species 

Ratio 

Species 

Ratio 

hc 2 n+ 

-4.5 

C 10 H 2 

- 1.2 

N 2 H+ 

-0.7 

C 4 H“ 

-0.4 

HC 8 N 

-3.1 

CH 3 C 3 N 

- 1.2 

c 5 h 2 

-0.7 

I-C 3 H 

-0.4 

hc 6 n 

-2.9 

H 3 C 4 NH+ 

- 1.2 

C 5 H 3 N+ 

-0.7 

c 6 - 

-0.4 

hc 4 n 

-2.9 

C S H- 

- 1.2 

Mg+ 

0.7 

Na + 

0.4 

ch 2 nh 

-2.4 

c 9 

- 1.1 

C 5 H- 

- 0.6 

hc 3 nh+ 

-0.4 

CCS 

- 2.2 

CgH 

- 1.1 

C 7 

- 0.6 

nh 3 

-0.4 

HNC 

-1.7 

c 9 - 

- 1.1 

c 7 - 

- 0.6 

ocs 

-0.4 

CioH 

-1.7 

c 9 n 

- 1.0 

n 2 

- 0.6 

CCN 

-0.4 

CioH- 

-1.7 

c 9 h 2 n+ 

- 1.0 

OCN 

- 0.6 

c 6 h 2 

-0.4 

hc 9 n 

- 1.6 

c 3 h+ 

- 1.0 

C S H 

- 0.6 

h 2 cn 

0.4 

nc 4 n 

- 1.6 

C 2 S+ 

-0.9 

c 7 h 2 n+ 

- 0.6 

CgH+ 

-0.4 

NC S N 

- 1.6 

c 7 h 

-0.9 

i-c 3 h 2 

0.6 

h 2 ccn 

-0.4 

nc 6 n 

- 1.6 

C 7 H' 

-0.9 

CNC+ 

0.6 

NH+ 

-0.4 

Cn 

- 1.6 

C 8 H 2 

-0.9 

CsN’ 

- 0.6 

c 6 

-0.4 

CH 3 C 7 N 

-1.5 

NH 

-0.9 

c-C 3 H 

- 0.6 

C S N 

-0.4 

hc 7 n 

-1.5 

HC 2 S+ 

-0.9 

nh 2 

- 0.6 

HNCCC 

-0.3 

H 3 C 7 N 4 - 

-1.4 

Cg 

-0.9 

HNO 

- 0.6 

c 6 h 6 

-0.3 

CH 3 C 5 N 

-1.4 

c 7 h+ 

- 0.8 

c 9 h+ 

-0.5 

C 5 H+ 

-0.3 

c 9 h 

-1.4 

C 10 H+ 

- 0.8 

ch 3 cn+ 

0.5 

CN' 

-0.3 

hc 3 n 

-1.3 

c 6 h 

- 0.8 

SiCgH 

-0.5 

HCNCC 

-0.3 

c 9 h 2 

-1.3 

ch 2 cch 

- 0.8 

CH 3 NH 2 

-0.5 

C 3 0+ 

mm* 

c 9 h- 

-1.3 

C 6 H- 

- 0.8 

CH 3 NH 3 

-0.5 

NS 

mm-: 

C 10 

-1.3 

Fe+ 

0.8 

c 4 h 

-0.5 



hc 5 n 

-1.3 

C 7 N 

-0.7 

NH 4 

-0.4 



cr 0 

-1.3 

Cg 

-0.7 

ch 3 c 4 h 

-0.4 




12 




















Table 5: Species with abundance ratios log(X new /X 0 id) smaller than -0.3 or larger than 0.3 at 10 6 yr, and 
with an abundance in at least one of the models larger than 10~ 12 compared with the total proton density. 


Species 

Ratio 

Species 

Ratio 

I-C3H2 

1.6 

hc 3 n 

0.5 

HC 8 N 

-1.5 

SiS 

-0.5 

SiCgH 

-1.4 

C 3 S 

0.5 

CCS 

-1.4 

N 2 0 

0.5 

hc 6 n 

-1.2 

ch 2 chc 2 h 

-0.4 

NO 

1.1 

c 3 0 

0.4 

nc 8 n 

-1.1 

Si0 2 

-0.4 

SiC 6 H 

-1.0 

NO+ 

0.4 

H 2 CCN 

-1.0 

OCN 

-0.4 

nc 6 n 

-0.9 

c 4 h 

-0.4 

CgHg 

-0.8 

h 2 nc+ 

0.3 

c 3 n 

0.8 

1-C 3 H+ 

0.3 

ocs 

-0.7 

nh 2 

-0.3 

CNC+ 

0.7 

HNO 

-0.3 

HNCCC 

0.7 

C 3 

0.3 

HCCNC 

0.6 

N 

0.3 

hc 4 n 

-0.6 

1-C 3 H 

-0.3 

c-C 3 H 

-0.6 

c 4 h 2 

-0.3 
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Fig. 2.— Mean logarithmic difference between the modeled and observed abundances over all observed 
species (see text) as a function of time using kida.uva.2011 (solid line) and kida.uva.2014 (dashed line). 
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Fig. 3.— Absolute logarithmic difference between modeled and observed abundances for each species at 
the best time of each model using plus signs for kida.uva.2011 and diamonds for kida.uva.2014. Horizontal 
lines show the differences between the old and new values. Red lines indicate that the difference has been 
increased by the new network while blue lines indicate that it has been reduced. Boldface species are those 
with abundances improved by the new network. 
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Fig. 4.— Abundances of CH 3 CHO and CCO as a function of time computed with the kida.uva.2011 (grey 
lines) and kida.uva.2014 (black lines) networks. The vertical grey and black lines correspond to the best 
times constrained with the old and new networks respectively. 
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